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Abstract 

An empirical fit is described to measurements of inclusive inelastic electron-deuteron cross sec- 
tions in the kinematic range of four-momentum transfer < < 10 GeV'^ and final state invariant 
mass 1.1 < W < 3.2 GeV. The deuteron fit relies on a fit of the ratio Rp of longitudinal to trans- 
verse cross sections for the proton, and the assumption Rp = Rn- The underlying fit parameters 
describe the average cross section for a free proton and a free neutron, with a plane-wave impulse 
approximation used to fit to the deuteron data. Additional fit parameters are used to fill in the 
dip between the quasi-elastic peak and the A(1232) resonance. The mean deviation of data from 
the fit is 3%, with less than 4% of the data points deviating from the fit by more than 10%. 

PACS numbers: 25.30.Fj,13.60.Hb, 14.20 Gk 



* Electronic address: |bosted@jlab.org| 
'''Electronic address: |christy@jlab.org 



1 



I. INTRODUCTION 



Empirical knowledge of the inclusive electron-deuteron cross section in the nucleon reso- 
nance region is important input for many research activities in nuclear and particle physics. 
The most important examples include calculations of radiative corrections to cross sections, 
extractions of spin structure functions from asymmetry measurements, determinations of 
structure function moments, and determinations of the vector coupling in models of low 
energy neutrino-nucleon cross sections. The latter is of particular importance since the 
quality of low energy neutrino-nucleon cross section models will become one of the largest 
uncertainties in the extraction of neutrino oscillation parameters from future long-baseline 
experiments. 

In this paper we will describe a fit to precision electron-deuteron inelastic cross sections 
in the resonance region for negative four-momentum transfer < < 10 GeV^ and final 
state invariant mass 1.1 < < 3.2 GeV. Among the advantages over previous fits [l, 2] are: 
inclusion of photoproduction and low data points; inclusion of several new experimental 
results; and the bulk of underlying fit is to the free nucleon (average of free proton and neu- 
tron), with Fermi motion consistently taken into account in a Plane Wave Impulse Approx- 
imation (PWIA). The latter provides for a parameter- free way of describing the broadening 
of resonance peaks with increasing three-momentum q. The average nucleon transverse cross 
section fit form is the same as the recent fit to electron-proton data of Ref . , which uses a 
set of threshold-dependent Breit-Wigner forms for all resonances. The principal difference 
is that the present fit is to the transverse portion of the cross section only, due to a lack of 
sufficient virtual photon polarization e range in the presently available deuteron data set. 
The assumption was made that the ratio of longitudinal to transverse cross sections, i?, is 
the same for the free proton and neutron, as supported by a recent analysis of data with 
W > 2 GeV {41, and also supported by the limited study of the present analysis. With 
the above assumptions, the "dip" region between the quasi-elastic peak and the A(1232) 
resonance is under-predicted at low Q^, possibly because Meson Exchange Currents (MEG) 
and Final State Interactions (FSI) have been ignored An additional empirical function 
was used to fill in the missing strength in the "dip" region. 
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The overall structure of the fit can be summarized as: 




where cr^{W, Q"^) is the transverse cross section for a deuteron, cJdipiW, Q^) is the dip region 
parametrization (6 free parameters), and the integral is over Fermi momentum k of the 
average free nucleon transverse cross section crJf{W, Q^) (36 free parameters). By simplifying 
the 3-dimensional integration to a 1-dimensional integration along the direction of the virtual 
photon, it was possible to simultaneously fit all 42 parameters in a single gradient search 
minimization. The minimization was done with respect to an ensemble of approximately 
15,000 data points from 9 experiments. 

The following section defines terms and kinematic variables. Section III describes the 
data sets used. Section IV gives details of the functional form used for a^ipiW, Q"^) and 
aJfiyV, Q"^), and how the Fermi-smearing integral was simphfied. The fit parameters are also 
fisted in this section. In Section V, we discuss various features of the results. 

II. DEFINITIONS AND KINEMATICS 

In terms of the incident electron energy, the scattered electron energy, E\ and the 
scattering angle, 9, the absolute value of the exchanged 4- momentum squared in electron- 
nucleon scattering is given by 



with Mp the proton mass, u — E — E', and the small terms involving the electron mass 
squared have been neglected. 

In the one-photon exchange approximation, the spin-independent cross section for inclu- 
sive electron-nucleon scattering can be expressed in terms of the photon helicity coupling 
as 



Q'^i-qr^iEE'sin'^, 



(1) 



and the mass of the undetected hadronic system is 



= Ml + 2MpU - Q\ 



(2) 



da 



V[al{W,Q^)+ea^^{W,Q^)\, 



(3) 



dndE' 
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where ajj (a^) is the cross section for photo-absorption of purely transverse (longitudinal) 
polarized photons, 

r = ^ ^— (4) 



is the flux of transverse virtual photons, and 

(5) 



l + 2(l + -)tan^- 



is the relative flux of longitudinal virtual photons. All the hadronic structure information 
is contained in aj^ and cr^, which are only dependent on W and Q^. We use the definition 
Rn = o'n/^'n- W^e use the subscripts p, n, A^, and D to refer to proton, neutron, average 
nucleon, and deuteron respectively, where the deuteron cross sections are defined to be per 
nucleon rather than per nucleus, following the high energy convention. We define = 

III. TREATMENT OF EXPERIMENTAL DATA USED IN FIT 
A. Description of the Data Sets 

The characteristics of the data sets used in the fit are summarized in Table HI The first 
reference (Ref. jsl) includes results from three early photoproduction experiments. We only 
used data from these experiments with beam energies above 800 MeV, because at lower 
energies the more recent photoproduction data of DAPHNE [3] has significantly smaller 
systematic errors. The JLab CLAS data fs] cover a wide kinematic range with many data 
points. The data are reported as values of the structure function F2, averaged over two 
different beam energies. Since the relative weight from the two beam energies (with differing 
values of e) was not given, there is a systematic error in the conversion to ctt that was taken 
into account in the total error bars. The early JLab Hall C data of Niculescu cover a 
similar kinematic range as the CLAS data, with less data points but smaller statistical and 
systematic errors. The more recent JLab Hall C data of EOO-116 [9] cover the high range 
with relatively few data points, which nonetheless have very small statistical and systematic 
errors (typically a few percent). To extend the low W region to even higher Q^, although 
with larger relative errors, we included the data of SLAC E133 [l^. To cover the lower 
region, we included preliminary data from two JLab Hall C experiments: E02-109* 13] and 



EOO-002 jl^. These experiments cover a wide range of e for each {W,Q'^) point with small 
statistical errors. Systematic errors were being finalized at the time this fit was done, so we 
used 2% as a conservative estimate in lieu of the ultimate errors for this experiment. 



Data Set 


Quin 


r)2 

^ Max 


^ Data Points 




(GeV2) 


(GeV2) 

V / 




Photoproduction (1972) [6] 








242 


Photoproduction (DAPHNE) \7} 





57 


CLAS [8j 


0.35 


5.9 


11725 


Early JLab [2] 


0.50 


4.2 


600 


JLab EOO-116 [9] 


3.6 


7.5 


288 


SLAC E133 [10] 


2.5 


10.0 


488 


JLab E02-109* [13] 


0.02 


2.0 


1435 


JLab EOO-002* [14] 


0.05 


1.5 


1445 


SLAC E140 [15] 


2.5 


10.0 


48 



TABLE L Data sets used in fit. The number of data points and the range are indicated for 
each data set. *The data from Refs. [l^. [l4| are preliminary. 



In order to constrain the fit at high W and Q^, where there are insufficient data from 



JLab, we included the DIS data from SLAC E140 [15|, and also added pseudo-data points 

n 

from the SMC [ll] fit to DIS ed structure functions. The pseudo-data were generated over 
the interval 2.4 < W < 3.2 GeV and 1.1 < < 10 GeV^. The errors used were those given 
by the SMC fit. 



The fit was found to be stable against the removal of any particular data set. 



B. Quasi-elastic subtraction 

Inelastic electron scattering on the deuteron can be divided into two distinct contribu- 
tions: quasi-elastic scattering (just proton and neutron in the final state), and inelastic 
scattering (one or more mesons in the final state). Since the goal of the present work is to 
fit inelastic scattering on the average free nucleon, we have subtracted the quasi-elastic con- 
tribution (if not already done by the experimenters) using the model described in Appendix 
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I. Data points for which the quasi-elastic fraction was greater than 30% of the total cross 
section were discarded, as well as all points for which W < 1.1 GeV. 

C. Longitudinal cross section subtraction 

The next step in the data treatment was to extract cr^ from each of the electroproduction 
cross section measurements. The was done using: 

(^D = (TD/il + eRD). 

As outlined in the introduction, we made the assumption that Rn = Rp, and evaluated Rd 
by Fermi-smearing both and crj from the proton fit of Ref. The Fermi-smearing 
procedure is described below. In practice, the Fermi- smearing had a very small effect for 
most {W, Q"^) points, so that, to a good approximation, Rd = Rp. 

IV. DESCRIPTION OF THE FIT 
A. Overview 

A gradient-search minimization code (MINUIT) was used to simultaneously determine 
all fit parameters by minimizing the value of x^? defined by: 

N 

1=1 

where the sum is overall experimental points with transverse inelastic cross section ai{Wi, Qf) 
and total statistical and systematic error 5ai{Wi, Q'j). To avoid a tedious iterative fit proce- 
dure, the Fermi-smearing integral in the model cross section a^iW^Q"^) was simplified to a 
single- dimension integral, as explained below. Starting values of the average nucleon cross 
section (jJfiW'Q"^) parameters were chosen to be the same as for the equivalent parameters 
in the proton fit of Ref. Some adjustments of the upper and lower limits on the parame- 
ters were needed to obtain the best ■ The functional form and starting parameters for the 
"dip" cross section adipiW, Q^) were obtained from detailed study of the fit residuals with 
the "dip" model absent. 
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In the next sections, we first describe the simphfied Fermi-smearing procedure, then the 
functional form of the free average nucleon cross section, along with the fit parameters 
obtained, and finish with the fit form and parameters for the "dip" region. 



B. Fermi-smearing 

The Fermi-motion of the nucleons in the deuteron was taken into account using a PWIA 



calculation and the Paris 



12j deuteron wave function 

a^iW, = j a^iW', {Q^)')^\k)Sk (6) 



where we made the approximations: 
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{W'f =\Md + v-\jM^ + k-^j - -{kf + 2q-k (7) 

and [Q'^y = Q^, where Md is the deuteron mass, M is the average nucleon mass, and u 
and q are the virtual photon energy and momentum, respectively. Since these equations 
basically boils down to the probability of finding a nucleon with longitudinal momentum k^, 
with the z axis chosen along q, we simplified the problem by determining 20 values of kl for 
which the integral over is close to 1/20, in the special case a]y{W', (Q^)') = 1. The 

corresponding average values of (fc^)j were also evaluated. The cross section is then 

20 

an{W,Q') = J2MW;,{Q'y)/20 (8) 

1=1 

where now 

{W[f = (m, + z/ - VM2 + (F).)' - (q)' - {k% + 2qkl, (9) 

This simplification is made possible by assuming that there are no off-shell cross section 
corrections, and that the k"^ terms are small enough that the integral over k"^ at fixed kz 
can be replaced by evaluating W at the average values of fc^. This approximation works 
well for ranges in W over which the cross section is slowly varying. This is the case for the 
kinematic region of the present fit, but is not advisable for smearing into regions where the 
elementary cross section is zero (i.e. W < M-|-M,r)- Hence, our fit is only valid for W > 1.1 
GeV. 

The numerical values of fc* and (fc^)j are listed in Table HTl 

Due to the rapid variation of the cross section near threshold, we used 200 bins instead 
of 20 for W <1.3 GeV. 
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i = 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 




0.0029 
0.0050 


0.0083 
0.0051 


0.0139 
0.0055 


0.0199 
0.0060 


0.0268 
0.0069 


0.0349 
0.0081 


0.0453 
0.0102 


0.0598 
0.0140 


0.0844 
0.0225 


0.1853 
0.0964 



TABLE II: The first ten values of kl and {k'^)i, in units of GeV and GeV^ respectively. The other 



10 values are given by kl'^^^ = and {k'^)i+iQ = {k"^) 



Free Nucleon Fit Form 



The fit form used to describe ajj (the transverse cross section for the average of a proton 
and a neutron) used the same functional form as Ref. 3|. The total cross section is defined to 
be the incoherent sum of contributions from resonance production (o"^) and a non-resonant 
background {a^^). The resonant cross section are described by threshold-dependent rela- 
tivistic Breit-Wigner shapes with Q^-dependent amplitudes for each resonance, such that 



a'^iW^, Q^) = J2 BW\W^) ■ Al{Q^ 



(10) 



The form used for the Breit-Wigner resonance shapes is given by 

■ptot-p7 



with 



K = 



{W^-Ml)/2M^, 
= {W^ - Ml)/2W. 



(11) 

(12) 
(13) 



Here, K and K'^'^ represent the equivalent photon energies in the lab and center of mass 
(CM) frames, respectively, while Ki and K^"^ represent the same quantities evaluated at the 
mass of the i*^ resonance, Mj. r*°* is the full decay width defined by 



-'tot 



(14) 



with /5j the branching fraction to the j decay mode for the i resonance and T* the partial 
width for this decay mode. The partial widths for single pion or eta decay were defined as 

L 

(15) 





2L+1 




.Pj \Mi_ 
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where the p^"^ are meson momenta in the center of mass, L is the angular momentum of the 
resonance, and Xi is a damping parameter. For two-pion decays, we used: 



The virtual-photon width was defined by: 





2 







cm\ \2 I y2 



cm\2 



+ x 



(16) 



(17) 



Since BW^{W'^) depends only on W"^, it was evaluated in 1 MeV bins in W from pion 
threshold to 5 GeV and stored in a look-up table for future reference. This significantly 
reduced the time needed for evaluation needed by the fitting code. 



i 


Mi 




U 










Ai{0) c\ 






1 


1.230 


0.136 


1 


0.145 


1.00 


0.00 


0.00 


8.122 5.19 


3.29 


1.870 


2 


1.530 


0.220 





0.215 


0.50 


0.00 


0.50 


6.110 -34.64 900.00 


1.717 


3 


1.506 


0.083 


2 


0.215 


0.65 


0.35 


0.00 


0.043 191.50 


0.22 


2.119 


4 


1.698 


0.096 


3 


0.215 


0.65 


0.35 


0.00 


2.088 -0.30 


0.20 


0.001 


5 


1.665 


0.109 





0.215 


0.40 


0.60 


0.00 


0.023 -0.46 


0.24 


1.204 


6 


1.433 


0.379 


1 


0.215 


0.65 


0.35 


0.00 


0.023 541.90 


0.22 


2.168 


7 


1.934 


0.380 


3 


0.215 


0.60 


0.40 


0.00 


3.319 





2.0 



TABLE III: Resonance parameters for states included in the fit. The branching ratios to single 
pion, double pion, and r] are denoted by /3i7r, (32n, and respectively. The assumed angular 
momentum is denoted by L*. Units of cross section are fib and all masses, momenta, and energies 
are in units of GeV. 

For the transition amplitudes the fit form utilized was 

The parameters for all the resonances are listed in Table [XTTl The variables Ai{0), c\, C2, 
and C3 were free parameters in the fit, while all other parameters were fixed to those used 
in the proton fit of Ref . ^ . 
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The non-resonant background was parametrized as 

= J2^'iCli6Wf'+'y^)/iQ^ + C'^)iC's+ciQ^+ciQ^) (19) 

i=l 

where 6W = W-Mp-M^, is the pion mass, and x' = 1 + {W^-{Mp + M^)^)/{Q^ + Cg), 
and the fit parameter Cq = 0.05 GeV^. The fit values for the other ten parameters are hsted 
in Table lYl 



i 


q 


CI 


CI 


CI 


CI 


1 


226.6 


0.0764 


1.4570 


0.1318 


-0.005596 


2 


-75.3 


0.1776 


1.6360 


0.1350 


0.005883 



TABLE IV: Non-resonant parameters as described in the text. Units of cross section are fib and 
all masses, momenta, and energies are in units of GeV. 

D. Dip region parametrization 

It was found that with the assumption that the PWIA "smearing" of free nucleon cross 
sections describes electron-deuteron scattering, there was always missing strength in the 
"dip" region between the quasi-elastic peak a.t W = Mp, and the low W side of the A(1232) 
resonance. This missing strength could be due to MEC and FSI in either quasi-elastic or 
inelastic scattering. Whatever the physical cause, we chose a purely empirical form that 
greatly improved the fit quality. This dip region additive correction term is given by: 

P = i.964z/0-298e[-(^-i-°^^)'/° ™] /(I + max(0.3, Q') /1.265)« (20) 

The unit-less structure function Fi is related to ar by 

Fi = {W^ - M^)aT /87T^ a{hcy (21) 

where a is the fine structure constant. 

V. FIT RESULTS 

The results of the deuteron fit are shown at seven representative values of as a function 
of W in Fig. [H along with the 1-a error band. In order to span a smaller range on the vertical 
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axis, we plot Fi rather than ctt- Note that the resonant structure, clearly visible at low Q^, 
is essentially gone for > 5 GeV^, due to the increasing influence of Fermi broadening and 
non-resonant background contributions. 

5.000 



1.000 
0.500 

0.100 

0.050 

0.010 
0.005 

1.25 1.50 1.75 2.00 2.25 2.50 
¥ (GeV) 

FIG. 1: Fit results for Fi/nucleon for the deuteron as a function of W at representative values 
of as indicated on the figure. Central values are shown as the thick lines, while the 1-a error 
bands are shown as the thin lines. 

Ratios of cross sections for each data point to the corresponding fit value are shown in 
Fig. [2] as a function of W for six bins in Q^. As shown in Fig. [3l 96% of the points lie 

11 




within 10% of the fit, 76% he within 5%, and 53% he within 3%. Overah, the agreement of 
data and fit is reasonably good at the 3% to 5% leveL It can be seen in Fig. [2] that there 
are two noticeable oscillations in the ratios at low W and Q^: this trend is also seen in the 
p.o.o„ m of Ref. S. Th. ,„a. ,„d,cate the need fo. add.iona, pa..a,„e.e. .o desc.be the 
photoproduction and very low data. 

In order to check the assumption that Rp = Rn, the data with e < 0.5 are plotted in 
gray (blue online), while those with e > 0.5 are plotted in black. No glaring discrepancy 
is seen between these two data sets. A more refined analysis will be performed once the 



preliminary data of JLab E02^109 [13| and EOO-002 [iJ] are finalized, and results from the 
very recent JLab E06-009 [161] experiment become available. The E02-109 and E06-009 
experiments were specifically designed to measure R^, and cover a large range of e at many 
specific values of (W,Q'^). 

n 

The ratio of the commonly-used fit of Niculescu [2|] to the present fit are also shown in 
Fig. [2l That fit tends to systematically lie above the data in the resonance region {W < 2 
GeV), and lies well below data at high W and low (outside the kinematic range of that 
fit). 

Since the present underlying fit is to the average free nucleon, we can use the proton fit 
to obtain predictions for the ratio of neutron to proton transverse cross sections (or equiva- 
lently the ratio F^/Ff), as illustrated in FigH^. Significant resonance structure is predicted, 
especially at low and in the region of the A (1232) resonance, for which the resonant con- 
tribution to F^/Ff is expected to be unity by isospin invariance. These predictions can be 
tested against the anticipated results of the JLab "BONUS" experiment which used 
tagging of low energy backward protons to "tag" spectator protons in electron-deuteron 
scattering (and hence isolate electron- neutron scattering). Predicted ratios of F^ / Ff that 
could be extracted from BONUS are shown in Fig|4)D at three representative values of Q^. 



VI. SUMMARY 



An empirical fit to inelastic electron-deuteron scattering has been performed which de- 
scribes available data reasonably well (3% to 5% level) in nearly all of the kinematic range 
than can be accessed at Jefferson Lab with up to 6 GeV electrons and photons: < < 10 
GeV^ and 1.1 < W < 3.2 GeV. The fit is useful in the evaluation of radiative corrections 
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1.2 
1.1 
1.0 
0.9 
0.8 
1.2 
1.1 
1.0 
0.9 
0.8 
1.2 

1-1 

X! 1.0 
O 0.9 
0.8 
1.2 



CO 



1.1 

1.0 
0.9 
0.8 
1.2 
1.1 
1.0 
0.9 
0.8 
1.2 
1.1 
1.0 
0.9 
0.8 



0.0<Q''< 0.2 
GeV^ 



H 1 1 — H 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1- 



i 1 1 1 1 1 1 1 1 HH 1 1 1 1 1 1 1 1 1 1- 




1^ 



0.2<Q''< 0.5 
GeV^ 



H 1 1 1- 



H 1 1 1- 



1^ 



0.5<Q''< 1.0 
GeV^ 



i 1 1 1 1 1 1 1 1 1 1 1 H^H 1 1 j 1 1 1 1 1- 



2.0<Q''< 5.0 
GeV^ 



1 — I — 1 1 I ii — — I — I — I — I — — I — I — I — I — — I — I — I — I — — I — I- 



GeV^ 



1.0 



1.5 



2.0 2.5 3.0 

W (GeV) 



3.5 



FIG. 2: Ratio of ed inelastic data to model as a function of W in six ranges of Q^. The gray (blue 
online) points correspond to e < 0.5, while the black points are for e > 0.5. Most of the points 
for W > 2.4 GeV and > i GeV^ are from the ed fit from SMC njl. The ratio of the fit of 
Niculescu [2j to the present fit is illustrated by the dashed curves. 

to experimental data, for extraction of spin structure functions from asymmetry measure- 
ments, and for the evaluation of structure function moments. Since the underlying fit is to 
an average nucleon, the results can be combined with a proton fit 3|] to obtain predictions 
for electron-neutron scattering in the resonance region. Suitably corrected for Fermi motion, 
these can in turn be used to make neutron excess corrections to nuclear structure functions. 

13], EOO-002 Q, and JLab E06-009 



Once the data from JLab E02-109 



16| are finalized. 



we plan to re-do the fit for ai and ctt separately, rather than making the assumption 
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400 



200 



T 1 1 r 



T 1 1 r 



T 1 1 r 



J I I L_ 



J I I L 



10 
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T 1 1 r 



15 



20 



FIG. 3: Frequency distribution for the deviations from unity of the ratios of data to fit. 



Rp = Rn- We also plan to use the results of the BONUS 17| experiment, which to first 



order will measure the ratio of electron-neutron to electron-deuteron scattering. 

FORTRAN computer code embodying the electron-deuteron fit described in this article 
is available by email request from the authors. The code includes the fit covariance matrix 
(or error matrix), and a subroutine to obtain the error on Ff. The code also includes the 

n 

proton fit of Ref. [3|, permitting the determination of electron- neutron cross sections from 
the deuteron and proton fits. The same code also includes a simple Fermi-smearing model of 
electron-nucleus cross sections for A > 2, using these deuteron and proton fits, as described 

n 

in Ref. [21|. The quasi-elastic model described in the Appendix is also included. 
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VII. APPENDIX I: QUASI-ELASTIC MODEL 

To model quasi-elastic scattering, we used the same PWIA Fermi-smearing prescription 
(based on the deuteron Paris wave function) as for inelastic scattering, except that the 
continuous inelastic cross section was replaced by a 5-function elastic cross aiW = M„. The 

n 

elastic cross section was calculated using the nucleon form factors of Bosted I18|, modified 

n n 

for off-shell effects using the prescription of Ref. [19|. Following Tsai [20], Pauli suppression 
was taken into account using the factor {3q/Akf)l — [{q/kf)"^]/ 12, for q < 2kf, and unity for 
q > 2kf, where q is the magnitude of q, and we used kf = 0.085 GeV. As shown in Fig. [5l 
the model works reasonably well near the quasi-elastic peak at two values of Q^. It can also 
be seen that the empirical "dip region" parametrization is useful to describe the region near 
W = 1.09 GeV at low Q^. 
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FIG. 4: Fit results for a) the ratio Ff /Ff and b) the ratio F^/F( as a function of W for = o.5 
GeV^ (soUd curves), = 1 GeV^ (dashed curves), and = 2 GeV^ (dot-dashed curves). F( is 
defined to be per nucleon. 
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FIG. 5: Comparison of the quasi-elastic model (dot-dashed curves), PWIA part of the inelastic 
model (dashed curves), "dip region" part of the inelastic model (short dashed curves), and their 
sum (solid curves) with F2 data from Ref. [8] at a) = 0.525 GeV^ and b) = 2.075 GeV^. 
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